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■ I. INTRODUCTION 

X) 

^ ] In recent years, a series of accurate calculations of N—d scattering observables below the deuteron breakup threshold 
(DBT) by using the Kohn variational principle (KVP) have been reported The resulting elastic n — d phase 

' shifts were found to be in strict agreement with those determined via accurate solutions of the Faddeev equations 
. (FE) j|] . The interest in this reaction has been revived consequently on the observation of large disagreements between 
^ ' theoretical predictions and experimental data for the p — d and n — d Ay and zTn polarization observables Since 
^ the solution of the N—d Schrocdinger equation can be obtained with high precision, the disagreements can be traced 
' back to the current models for the nucleon-nucleon (NN) and three-nucleon (3N) interactions. Therefore, the detailed 
T^J- , studies of this process can be used to get information on the microscopic nuclear interaction. 

Above the DBT, the FE have been accurately solved for the n — d reaction (see, for example, Ref. |^ and references 
C^l therein). For reactions where there are more than one charged particles in the final state, the proper inclusion of 
the Coulomb interaction is a difficult problem. On the other hand, the proton induced breakup reaction is more 
easily studied experimentally, and in fact, there is a large amount of high precision p — d scattering data available. It 
should be noticed that the effect of the Coulomb potential can be rather sizeable for various polarization observables, 
' in particular at small scattering angles. It is therefore of relevance to obtain accurate theoretical estimates of the 
I , observables of this reaction. 

Q ■ At present, different approaches are applied to calculating the p — d elastic observables above the DBT. The 
^ ' method used in Refs. consists of separating the long-range from the short-range Coulomb effects. Then, the 

^ ' FE in momentum space are solved by taking into account only the short-range part of the Coulomb potential, whilst 
I the long-range contribution is included by "renormalizing" the solutions found. This method has been pursued for 
• applications in both atomic systems and p — d scattering at i? = 10 MeV j8|J|] . In other approaches, the FE are solved 

using the Coulomb-Sturmian-space representation JlOt , or by directly solving the FE in configuration space . 

The KVP was discussed in Refs. p^-[l^ to treat breakup reactions without considering the case of two or more 
charged particles in the final state. In the past it has been used mainly for studying atomic processes (in particular, 
the scattering of electrons or positrons by hydrogen atoms). Above the ionization threshold the direct application 
of the KVP has however encountered various problems related to the infinite number of two-body bound states 
between particles of opposite charge and to the complexity of the boundary conditions to be satisfied [|l^,^ . Various 
alternative formulations have then been studied and the associated difficulties analyzed p7|-p^ . 

In this paper we reconsider the direct formulation of the KVP for calculating the elastic part of the S-matrix 
of a p — d collision above the DBT. In such a process there is no problem relating to two-body excited bound 
states. The long-range Coulomb interaction between the two protons introduces a distortion in the description of 
the asymptotic three-nucleon outgoing wave. We have proved the validity of such a principle, considering trial wave 
functions which satisfy "simplified" boundary conditions allowing for practical applications, even when using realistic 
nuclear interactions. The results for several elastic observables have been already reported in refs. j2n,El[. 



o 



The paper is organized as follows. Some definitions used later in the paper are given in Sect. O. In Sect. HI, a 



detailed description of the asymptotic behavior of the p — d scattering wave function (w.f.) is reported, while the proof 



of the KVP is given in Sect. IV. The use of such a principle in conjunction with the pair-correlated hyperspherical 



harmonic expansion technique is discussed in Sect. M, while the conclusions are given in the last section 
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II. DEFINITIONS 



Let us introduce the Jacobi vectors for three equal mass particles, given by 

[Af rj+rk\ 

Xi = rj-rk , yi = J-[ri — , (1) 



where i, j, A; is a cyclic permutation of 1, 2 and 3, and denotes the position of particle i. The distance between the 
pair jk and the particle i is 

di = KYi , K = J- . (2) 
For the p — d system, the non-relativistic total energy operator is written as 

^ = -^E^'+E ^''''('^0 + ^'"^ + , (3) 

i=l ]<k 

with V'^^(x,;) being the pair potential between the particles j, k and W^^ a 3N interaction term. Both V^^ and W^^ 
are short-ranged. In particular, V^'^(xi) and W^^ become vanishingly small when Xi > Qn, where ajv is a distance 
somewhat larger than the range of the nuclear interaction. Moreover, Vc is the Coulomb interaction operator (we 
disregard here other electromagnetic interaction terms), namely 

j<k ^ ^ 

where t„i.z/'2 is the z component of the isospin operator of particle m. 
For a three-body system, the hyperradius is defined as 

p = y/^i+Vi ' (5) 

and is independent of the permutation i = 1,2,3. The hyper angular variables are 

= {xi,yi,Oi) , (6) 

where the hyperangle 9i is defined by the relation 

cos.6i = Xi/p. (7) 
To describe a three-nucleon state with total angular momentum J, Jz, the following functions are introduced 

where the variables ,s [t) specify spin (isospin) states. The quantities Sa (T) and (f*) arc the spin (isospin) of the 
pair J, k and of the third particle i, respectively, and they are coupled to give Sq, (T). The index a stands for the set 
of quantum numbers 

a= {ia,La,Aa,Sa,^a} ■ (9) 

The indices {a, T, T} distinguish the various angular-spin-isospin "channels" . A generic three-body w.f. can be 
written as 

3 

M/ = ^V'(i), VW = V'(xi,yi), (10) 

i=l 

where each amplitude ip(i) corresponds to the i-rearrangement configuration of the three particles and it can be 
decomposed in channels as 
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(11) 



:,T,T 



The antisymmetry of ^ can be simply enforced by requiring that the amphtude -tjj^i) change sign under the exchange 
of particles j and k. This can be easily done by including only those channels having odd values of £a + Sa + T in 
the sum of Eq. (|l]) 

In the following, the projection on states with a definite value of the third component of the isospin of the three 
particles will be useful. Let us define 

E,^\tl = -l/2,ti = l/2,t'; = l/2) , (12) 
as the isospin state where particle i is the neutron. Therefore, each amplitude can be decomposed as 

3 



(13) 



i/=i 



where ipv{i) is the component of the amplitude ip{i) and particle v is the neutron. To expand the components V'i'(*) 
in spherical waves, it is convenient to use 

which are a complete set of angular-spin functions. The functions 3^a,T,r(i) and ya{i) are related to each other by 



1 „ 1 „ 1 - 



ya.jM') = y»i^) 



1 „ 1 - 

2 „ 1 „ 



(15) 



The components ip^ii) can be then expanded as the sum 



1 „ 



(16) 



where Fa^i, are radial amplitudes. Using Eqs. (p^), it is not difficult to find the relation between the amplitude -Fa,T,T 

defined in Eq. ( [TTI ) and the functions Fa^i,. 

The N — d system in the asymptotic region, where a nucleon and the deuteron are far apart, is described by the 
following "surface" function. 
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(17) 



where (jjd (^d) is the spatial-spin (isospin) part of the deuteron w.f. The quantum number L is the relative angular 
momentum between the lone nucleon and the deuteron, S is the spin obtained by coupling the spin j = 1 of the 
deuteron with the spin 1/2 of the third nucleon. The deuteron wave function ipd = 4>d{^i) x satisfies the following 
equation 



m 



(18) 
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III. THE ASYMPTOTIC BEHAVIOR 



The asymptotic behavior of the w.f. describing N — d scattering has been discussed by many authors, in particular 
for the n — d case [ p2|~p6| . Th e case of the scattering of three particles, with at least two of them charged, has been 
considered in Refs. 



d scattering 



]15| , |27| -|29|. In this section, we rederive the asymptotic behavior of the w.f. for a p 
process, in order to clarify some properties used to discuss the KVP in the following section. 

The asymptotic behavior of the amplitudes given in Eq. is analyzed by considering them as solutions of 
the three Faddeev-Noble equations (FNE) Q 



+ Vc + V''''{^,)-E 



(19) 



with i,jjfc cyclic. The 3N interaction W^^ has been disregarded since we are interested in the asymptotic solutions 
of Eq. ([19|), where at least one particle is very far from the other two. Let us consider the volume V{R) of the 
six-dimensional space where p > R {R ^ oo). V{R) can be divided into the following regions, corresponding to 
different ranges of values of the hyperangular variables fi. 

• 1+1+1 (breakup) region. Here all the particles are well separated from each other. This breakup region is 
denoted as Vb{R) = {p> R; Xi — + oo for i = 1, 2, 3}. 

• 2+1 regions. The part of V(i?) where the particles j and k are close enough to interact through the NN potential, 
while particle i is very far from them, is hereafter denoted by Vi{R). There are three such regions, corresponding 

to i = 1, 2, 3, namely Vi{R) = {p> R;Xi < flAr}. In Vi(i?), we have yt > R but the distance Xi between particles 
j and k is of the order or lesser than the range cat of the nuclear potential. As R increases, the range of the 
hyperangle 9i subtended by the region Vi{R) decreases, and in the limit i? — > oo it reduces virtually only to the 
value 9i — n/2. 

• Transitions regions. The regions between V&(i?) and Vi{R) are denoted as Vt,i{R), i = 1,2,3. They can be 
defined as Vt,i{R) = {p> R',Xi < Tp^^-.Xi > oat}, where T and ^ are suitable constants {0 < p < 1, as wiU be 
shown below). In this region 



cos 6i 



P 



< Tp' 



, for p — > oo 



(20) 



Therefore, for p — s- oo the hyperangle extension of this region reduces again to the value 9i — tt/2. Such regions 
were considered for the first time by Merkuriev ||lj,|2^ and should allow for a smooth transition between the 
solutions in the regions Vb{R) and Vi{R). 

Let us study the asymptotic behavior of the solution of the FNE given in Eq. (^9|) in the above described regions. 

A. The Vb{R) region 

In the breakup region Vb{R) the nuclear interaction can be safely disregarded and the FNE for ip{i) becomes 



m ^ " 



Vc~E 



ijj{i) = . 



(21) 



By rewriting this equation using the p, fit coordinates and disregarding terms vanishing faster than 0{l/p) as p ^ oo, 
the solution reduces to ETIl 



Aa.T.T{Oi)ya.T.T{i) 



+ 



g-iQp+ixln2Qp ^ 



Ba,T,r{di)ya,T,r{i) , region Vb{R) , 



(22) 



4 



where E = ti^Q'^/m and 

m 1 + Tj 2 1 + Tfe 2 

*"2Fq^^^" 2^' ^^^^ 
^ I— 1 

is a dimensionless operator. The functions Aa^T,TiOi) and i?a,T,T(^i) are the so-called breakup amplitudes and are 
fixed by the dynamics of the process. For example, the cross section for the process p + d^p + p + n is proportional 
to lA^.T.rm''. 

Alternatively, 'ip{i) can be decomposed as in Eq. ([l^), where the neutron is explicitly labelled. The corresponding 
components a-^e 



M^) = :i ( E^k^(^03^o»j 

Y.^l.my"(')] ' region V6(i?) , (24) 



g-iQp+i?7(e„)ln2Qp 



where ri{9y) represents the isospin matrix elements of the operator x a-nd is defined as 

rjyp 

2ft QcosO„ 

The breakup amplitudes depending on the neutron label v are related to the previous ones by the relations: 



Al,m = -^A^^.^.ie.) - -La„_i,i(0,) + -La„,o^. (0,) , 

and analogously for i?^ j. The superscript i of the functions A^^ ^ recalls that these amplitudes are the appropriate 
combinations of functions A^ji'^t entering the expression (|2j) of ?Ai/(i). 

B. The Vt,i(-R) region 

In the following analysis we will fix our attention to the term describing outgoing breakup waves, i.e. proportional 
to exp(i(5p). When cos0,y 0, namely going into the regions Vt,j/(-R), the term exp(— i?7(0i,) In 2Qp) in Eq. ( [2^ ) 
oscillates increasingly. Moreover, it is no longer appropriate to disregard terms like l/(pcos0i/)^ in the construction 
of the solution. In this region the FNE satisfied by ■;/'(*) is still given by Eq. (pi]), since Xi > aj^. We can distinguish 
the following cases. 

1. Component ?/'«(*) (particle i is a neutron). The FNE satisfied by this component is 



vi. + v.? I + - - ^ 



Mi) = . (27) 



Since we are in a region where Xi <^ yi ~ p, the general solution of the above equation can be written as 



V'i(«) = X! / d,LpCa,i{^p)—Fi^ 



,qxA-h+(py{)y^(i). (28) 



where 

q = QcosLp, p = Qsiiiip, (29) 



5 



Ft[rj^qx) is the regular Coulomb function and 

hliz) = z(-yL{z) + -ijL{z)) , h+{z) ^ e'(--W2) for z ^ oo , (30) 

with ji and yL spherical Bessel functions. The function Ca,i{<f) entering Eq. (p8|), whose dependence on Q is implicit, 
is a smooth weight function, strictly related to the breakup amplitude Aa^T.riSi)- It goes to zero for (ys ^ or tt/2. 
In order to analyze the asymptotic form of the solution it is convenient to decompose the Coulomb function as 



Fi{ri, z) = Ti{'q, z) sm[z + /3<>(ry, z)^ , 

where J^i{ri, z) and Pi{ri, z) have the following asymptotic behaviors for z ^ max(l, ry^), 

Ti{r],z) = 1 , 

/3^(7?,z) = -r/ln22-£7r/2 + cr^(g) , 
ai{q) being the Coulomb phase-shift. For z <^ min(l, l/rj) we have instead 



2^e-'^"/^|r(^+l + ir;)| 
T{2i + 2) 



(31) 



(32) 
(33) 



(34) 
(35) 



where F is the gamma function. The decomposition (^) is somewhat arbitrary. However, we require that both J- 
and f3 always be finite. Thus, z„ + Pe^i], Zn) = mr for n = 0, 1, . . ., where z„ are the zeros of the function F((r], z). An 
example of the decomposition is shown in Fig. (VI) for a few values of 77 and i. 

In the large p limit the integral in Eq. ( p^ ) can be evaluated by the saddle-point approximation. Introducing the 
above decomposition for the Coulomb function and the relations Xi — p cos 9i and yt — p sin 9i , the leading term, 
given by the saddle-point ip = 9i, is 



27r 



,37ri/4_ 



c 



2i cos 6*, sin 6*, 



?/(6'i),(5pcos 



X exp 



COS^ 6i ] — iLa — 



Y.^^^ [v{0^).Q peas'" e, 

iQp-iT^{e,)\n{2Qp) 



X exp 



i^Q- - icr<;„(Qcos( 



1 



x^L,»(^«) — exp 
p2 



3^a(i) , 



(36) 



where in the last passage we have defined 
A] 



2i cos 9j sin 9j 



exp( -177(6',;) In (cos 9i) - i(£a + icr^ ^ (Q cos 6*., 



(37) 



The next term of the saddle-point expansion is ~ p In the region Vt,i{R) we have cos6'i ^ so that rf'{9i) ^ 1. 
Therefore, if 



Qpcos"9i > if[9i) , 



or equivalently 



Xi » Pc 



(38) 



(39) 
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the functions T and /3 reach their asymptotic values given in Eqs. ( |32| ) and (|33[). Then, '^iii) reduces to the form of 
Eq. ( p^ ) with j given by Eq. (|37|). When Eq. ( |39| ) is not satisfied, corresponding to the situation in which the two 
protons approach each other, the functions T and /3 start deviating from their asymptotic behavior and the breakup 
w.f. cannot be cast into the form ( p4[ ) and (|2^). The extreme situation corresponds to the condition 

1 



or equivalently. 



2r 



gp cos^ 61, < , (40) 

< = ^ w 60 fm , (41) 

and the functions T and /3 reach the behavior given in Eqs. ( |3^ ) and (p5[). Thus, — > e~'^™*^ Qcosdi _^ q gjjj(.g 
Xi/p — cos 9i 0. Also the successive terms in the saddle-point expansion are zero. This corresponds to the intuitive 
fact that the probability of finding two protons (in this case, labelled as particles j and k) at any finite distance 
Xi — r.jk is negligible in the asymptotic region R — > oo. We arrive at the important result that the amplitude V'j(*) 
goes smoothly from the behavior given in Eq. (^) to a vanishing value when going from the region where Xi S> Pc to 
the one where Xi <^ p'^. 

The same conclusion holds if the regular Coulomb function Fi{r},z) is replaced by the (regular) solution of the 
two-body Schroedinger equation 



2r] £{£+!) mV^^ix^) 



F/{ri,q,Xi) = , z = qxi 



(42) 



where, for the matter of simplicity, we have considered here a central NN potential V^^{xi). Outside the range of 
we have F/ Fe + tan Se Gi where G is the irregular Coulomb function and 5i the phase shift induced by the 
NN potential. However, for q — > or 77 — s- oo it can be shown that G — > exp(7rry) and tanSi — > exp(— 27r?7). Therefore, 
when Xi ^ p'^ again F^ — > 0. 

2. Component V'j(*) (the neutron is labelled as j). The equation satisfied by '4'j{i) in this region is 



m 



E 



(43) 



In Vt,i(i?), Xj ~ Ki/i and therefore the solution of the above equation can be written in the form 



dip Ca.j{(p)—je^ (qxi) 



1 r 



G, 



py^) + lFLJ 



3^a(*) . 



(44) 



where Gl is the irregular Coulomb function and ji is a regular Riccati-Bessel function. As has been done before, this 
latter function can be decomposed as. 



ii{z) = Ji(z) sin(z + Ci{z)) , 
where Ji(z) and C,i{z) have the following asymptotic behavior for z 3> 1, 



and for z <C 1, 



Ji{z) 



Mz) 



1 , 



-£7r/2 , 



(2^ + 1)!! 



Ci{z) = . 

Note that Ji{z) = J^^(0,z) and Qiz) = /3f(0,z). Moreover, Jo(^) = 1 and Co(z) = 0. 



(45) 



(46) 
(47) 



(48) 
(49) 
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Following a procedure similar to the one adopted in point 1), the following expression for the component ^Jj{i) in 
Vt.iiR) is found, 



ipj (i) = C{p) ^ Ji^ [Q cos e^Xi) exp 



where 



In Eq. (50), we have defined 



Therefore, for 



Kia. (Q COS 9tXi 



„3iri/4 



„ . ^ exp -i{ia + La) - +iae^{Q8ine,) 

2i cos Oj sm Oj \ 2 



1 / UIG 

C{p) ^ — exp iQp - i— 27— In 2Qp 
p2 \ zn LjK 



(50) 



(51) 



(52) 



> Pn = I ^ 



(53) 



where the behavior given in Eqs. (^) and ( ^7|) for Ji{qxi) and Q{qxi) is valid, the component V'j(*) coincides with 
Eq. ( p^ ) [note that cos9j ^ k in Vt,i{R)]- 

If a;i <C Pn, the amplitudes of the channels having > become vanishing small since Je^ oc (x»/p„)2^". This 
corresponds to the intuitive fact that the probability of finding two nucleons (in this case, labelled as particles j and 
k) with a vanishing relative velocity and relative orbital angular momentum greater than zero is negligible in the 
asymptotic region 00. Therefore, 



M^)^C{p)Y,6,^,oAl,{l)y^{i) 



X., < p„ 



(54) 



3. Component ipk{i) (the neutron is labelled as k). The solution for this case coincides with that discussed in point 2. 

Finally, we can estimate of the value of the constant p defining the region Vt^i{R). In the n — d case, the appropriate 
value is /i « 1 /2 [ see Eq. (|5^)]. In the p — d case, a value of /i « 3/4 can be chosen in accordance with Eq. (p9|). As 
shown in Rcf. [|lj[ , the presence of those regions does not alter the validity of the KVP for the n — d case. Therefore, 
it can be expected that such regions do not play any important role also in the p ~ d case, as will become evident in 
the next section. 



C. The V^{R) region 



In this region the nuclear interaction between particles j and k cannot be disregarded, whereas particle i is far from 
the other two. The solution consists of an elastic p — d amplitude ip'^^\i) plus a breakup amplitude tp^''\i), namely 
■0(i) = ?/;('^)(i) + The amplitude ip'^'^^ has the form 



LS 



aLs'Hti'no,qodi) + bLs'HL im^Qodi) 



where ^^si^) given in Eq. (|lj), the distance di in Eq. (||) and 

GLiva^qadi) ± ii^L(?7o, go^i) 



qadi 



The functions F and G are the regular and irregular Coulomb functions and 

Vo 



(55) 



(56) 



(57) 



is the usual Coulomb parameter, with the wave number go satisfying the relation 
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-B, + j^—ql=E. (58) 
(k)^ m 

The coefScients qls {^ls) a-re the ampUtudes of the corresponding incoming (outgoing) waves. In fact, Ti,^ w 
exp(±igo'ii + . . .) in the hmit of large qodi values. 

The amplitude ip^^"^{i), when particles j and k are protons, again vanishes for p — > oo due to the Coulomb repulsion, 

namely '>pf'\i) = 0. When one of the particles is a neutron, it takes the form 

V'(''H*)=5]<5,„,o0°(a:.)^L,.(|)C(p)5^o(z) , ,. = jovk, (59) 

a 

where is the s-wave solution of a inhomogeneous zero energy two-body Schroedinger-like equation |^,^ derived 
from Eq. (O) by taking into account the coupling between 7/1(1) and ■(/'(j) + i){k) induced by the NN potential. The 

amplitude V^^^j, j,(i) vanishes for £a > due to the centrifugal repulsion. For the present discussion of the KVP the 
precise form of the functions cj)'^ is not important. They are normalized so that 1 outside Vi{R), namely when 

Xi ^ a^r, in order to match the result obtained in the transition region. The derivation of the equation satisfied by 
(j)'^ is reported in the Appendix. 



D. The Vj{R)+Vt.jiR) regions 



In these regions yj — > 00 and Xj/yj — > and therefore we can approximate Xi ~ Kyj and yi 



component ipiii) satisfies again the uncoupled FNE given in Eq. (|21 
Eq. (§) 



Moreover, cos^i '--^ k or 



Mi) -> Cip) J2 AU^)yc.{i) , region V, (i?) 



yj/2. Then, the 
Tr/d) and from 



(60) 



In fact, here (for p 00) we always have Qpcos^Oi ^ rf'{6i). The same conclusion is reached for V'fc(i)- Both 
components can be obtained directly by taking the limit Oi — > 7r/6 in Eq. (p^. 

The component '0j(*) is more difficult to obtain. In fact, since here Xj can become very small, we have to fully 
take into account the distortion due to the Coulomb potential. This component can be calculated by expressing it in 
terms of the coordinates Xj and yj and by rewriting Eq. ( ^l|) as 



m 



e 



E 



(61) 



Now '0j(i;xj,yj) can be determined by following the same procedure as in Sect. IIIB. In particular, it has been 
shown there that for Xj <^ p'^, V'jX*) vanishes due to the Coulomb repulsion between the protons. 

The behavior of ip^i) in Vk{R) and Vt,k{R) can be analyzed in the same manner. Here the components ipi{i) and 
ipjii) are given by Eq. ( |60| ) and can be obtained directly as the limit of the solution valid in Vb{R). The component 
ipkii) will be distorted by the Coulomb repulsion and will vanish inside Vk{R)- 



E. Summary 

It is straightforward to show that the behavior of the incoming wave oc exp(— iQp + . . .) has similar properties. 
Therefore, the behavior of the amplitudes '>p{i) in the asymptotic regions can be put in the form. 



LS 



,iQp-ixln(2Qp) 



Aa.,T.,T{p,^i)ya,T,T{i) 



-iQp+iX ln2Qp 



Ba.T.T{P-,^i)ya.T.T{i) 



(62) 



,a,T,T 
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The "distorted" breakup amplitudes Aa,T,T and Ba,T,T have been written as functions of (p, fii) taking into ac- 
count the modifications induced by the Coulomb, centrifugal and NN poten tials. In Vb{R) ^ clearly Aa,T,T{pT^i) 
Aa,T,T{di) and Ba.T.riPj^i) — > Ba,T,T{6i)- Moreover, as discussed in Sect. IIIBj and [II C , in the region Vi{R) and 
Vt,i(i?), a good approximation for the (isospin-projected components of) Aa.T.T{p,^i) should be 



lL,fc(p,f^.)-.^l"^(a:.,e.)lL,fe(^^.) , 
where the following functions have been introduced, 

T^^''J{xi,9i) = Ti^ {jiiOi), Qpcos^ exp if3e^ (^T](ei),Qp 



(63) 
(64) 
(65) 



) cos^ 6, 



'ir]{0i)lii{2Qpcos^ Oi) +i^Q- - icr^^ (Q cos 6ii) 



.f1"^(x„0O = (5,^,o{<lt{x^) - 1) + l) Ji^iQpcos'e,^ 



X exp 



) cos^ 9i 



TT 



(66) 



(67) 



In the regions Vi{R) +Vt,i{R) the divergent behavior of x for cos^i — > is completely canceled by terms in Aa.T.T- It 
can be checked that outside of the region Vi{R) + Vt,i{R) the functions A reduce to the breakup amplitude A, since 
in Vb{R) J-"*-^-* (a;i, 6*^) — > 1 (z = c, rt). The factor J^^'^^ takes into account the effect of the Coulomb and centrifugal 
potential in Vt,i{R) and goes to zero in Vi{R). Conversely the factor JF(") takes into account the effects of the NN 
potential and the centrifugal barrier. Since for £a > Ji^ vanishes in Vi(i?), the influence of the NN potential has 
been taken into account only in the la = Q wave. 



IV. THE KOHN VARIATIONAL PRINCIPLE 



The process we are interested in is a collision between a proton and a deuteron, namely, the reaction 

(p + d) -> (p + d) + (n + p + p) . (68) 



If the incident state of the p — d system has quantum numbers Lq, and J, the elastic amplitude in Eq. ( |62[ ) is 
characterized by 



(69) 



The coefficient aj^s is usually renamed as —SLgSo,LS: ^LoSo,ls being the elements of the elastic S-matrix. Since the 
breakup part should not include incoming waves, we must have 



Ba,T,T{Xi,9i) = 



(70) 



The Eqs. (H-Jfl) give the boundary conditions to be satisfied by the w.f. describing the process (|68|). Explicitly, the 
asymptotic amplitude is 



LS 



JQp-ixln{2Qp) 



Aa,,T.,TiP,^^)yo.,T,Tii)] , P>R 



'■T.T 



Let us now consider the integral 



(-) 



(71) 



(72) 
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where Qn stands for the integration over cPxi cPyi of the six-dimensional volume with p < R {R ^ oo) and vl/(~^ is 
the time-reversed w.f. In the end the limit i? ^ oo is taken. The function 5* is a trial w.f. which will be chosen to 
satisfy the boundary conditions. Without loss of generality, the trial w.f. can be decomposed into three amplitudes, 
^' — V'C*)- The asymptotic behavior of tp{i) can be written as 



LS 



iQp-ixlog 2Qp 



J2 ^°''T,T{p,^i)ya,T,Tii)\ , P>R- 



(73) 



> Q,T,T 



The overlined quantities are free parameters which represent approximations to the exact values. It should be noticed 
that an approximate expression x of the operator x has been introduced in the above equation. 

The important point is that the validity of the KVP for the elastic part of the S-matrix can be proved also in 
this case. However, the "choice" of x is not completely arbitrary, since the matrix element — E)\'i') has to be 

finite. Moreover, in a practical calculation the convergence of the elastic S-matrix should be checked whenever an 
improvement is made to 

For example, in the case considered in Sect. ^ the operator x is represented as an expansion over the Hyperspherical 
Harmonic functions and therefore has not any singular behavior for 9i i = 1,2,3. Moreover, x ^ X the 

number of terms included in the expansion basis is increased. Such a procedure may result in a slow convergence of 
A towards the exact A, in particular for 9i 7r/2, i — 1,2,3. However, this problem has no practical consequences 
for the calculation of the elastic S-matrix, as will be shown in Sect. 

Contributions to the integral / come only from the differential operators present in H. By integrating by parts, / 
takes contribution only from the hypersurface [O] a,t p = R. Let us denote with [ilj,] the parts of [Q] corresponding 
to a particular asymptotic region Vx{R)- Therefore, 



h 



(74) 



where 7^, It.i and lb are the contributions coming from [fii], [flt.i] and [fib], respectively. On the hy persurface [Q], the 
trial and exact wave functions have reached the asymptotic behavior given in Eqs. (|7l|) and {\[^. In [n,], only the 
elastic part contributes, with the result 



T 2i n'f 

qa m\ 



(75) 



In fact, the breakup part goes like p^^/^ and gives a contribution cx R~^/^ which vanishes in the limit R — > oo. The 
contribution It^i coming from the region \^t,i\ is also zero. First of all, only breakup waves can contribute since the 
elastic wave vanishes outside Vi{R). Moreover, the "dimension" of the surface [^t.i] is proportional to whereas 
the integrand is proportional to 1/R^. Thus, for i? ^ oo the integral in [f2t,i] behaves oc 0. Finally, in the 

case in which only outgoing breakup waves are present in both the trial and exact w.f. the contributions of the two 
terms in Eq. ( [T^ ) cancel each other. 
In [fib], lb reduces to 



//, = / dfl 

[rib] 



*ioSo ^ioSo - ^LoSo^P^ *LoSo 



(76) 



and, using the asymptotic behaviors given in Eqs. (|62|) and (|73|), lb oc 1/R for i? ^ oo. 

Summing up all the contributions and remembering that (H — i?)'I'i(,,5„ = the exact expression is obtained: 



^LqSq.LqSo ^ LqSq,LqSo 



mqQ_,i-) 



{H-E)\eLoSo) , 



(77) 



where eloSo — ^ioSo ~ ^LoSo is the "error" introduced in the trial w.f. It can be noticed that Eq. ( [ttI) is variational 
in character, since the last term is quadratic in the error function eloSo- Therefore, the functional 



/■:t7(-) 



LqSq 



{H-E)\^iJLo 



Sol 



(78) 
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represents a second order estimate for the exact diagonal element SloSo,LoSo- The search of stationary solutions 
through the equation 

S[SloSo,LoSo] = > (79) 

provides the optimum estimate. Here S denotes the variation with respect to all the trial parameters entering '^LoSo, 
including S and A. Moreover, the KVP can be readily generalized to obtain variational relations for the non-diagonal 
matrix elements SlxSi,L2S2 of the elastic collisional matrix. 

Clearly, the principle is meaningful only if the two matrix elements entering Eq. ([77^) a re both finite. This is not 
evident since for p > R both and e have oscillatory behaviors. However, since Eq. (|77|) is exact and S and S are 
finite, it is sufficient to check that one of the matrix elements is finite, in particular {^\{H — This imposes 

some constraints on x, which usually has to be chosen coherently with the expansion basis used to construct ^' (see 
Sect. ^ for a practical example). 

Applications of the KVP usually start with the expansion of the w.f. in a certain basis which is truncated at some 
level. The stability of the functional ( |7^ ) is studied in terms of the dimension of the basis and in terms of the nonlinear 
parameters (if any). Spurious resonances or spurious solutions in the convergence procedure could appear due to a 
zero eigenvalue present in the operator H — E Q . These spurious solutions have been studied for a long time and it 
was shown that the KVP in its complex form, as presented here, is well-behaved |^^. Spurious solutions could appear 
in this case too, but these occur for unrealistic values of the nonlinear parameters, which are in general present in the 
form of a regularizing factor or inside some exponent Js^ . Below the DBT, a practical rule to verify the presence of 
spuriousities is to apply the principle to the reactance matrix K and its inverse (inverse KVP) and verify the 
relation K K^^ = I. Similarly, for energies above the DBT, one could apply the variation principle for the T-matrix 
(which can be derived by following an analogous procedure to the one described here) . With two independent estimates 
of (the elastic part of) the S~ and T-matrices one can test, for example, the relation S = I — 2iT. Besides these 
anomalies, which constitute an inherent problem of the method, the trial w.f. must be flexible enough to approach 
the exact one in the asymptotic and the internal regions. In other words, the error introduced due to the approximate 
representation of the operator x in the truncated basis should be reduced accordingly with the improvement in the 
description of the internal part of the wave function. 

The situation is rather similar to that which occurs when an approximate form 0^ of the deuteron w.f. is used 
in Eq. ( p^. It is easy to show that the KVP is still valid, but the "error" term also includes a quadratic term in 

~ 4'd P4|]- In general, what happens is the following: when the inaccuracies of or x are "small", the resulting 
errors in the elastic S-matrix elements are of second order, and can usually be disregarded. If the inaccuracies are not 
small enough, the expansion basis {Q} used to describe the internal part of the w.f. will spread out in the asymptotic 
region in an attempt to remedy the defects of the and x. If this happens, the convergence with respect to {Q} 
may be extremely slow, or it could even happen that no convergence could be found. 



V. A PRACTICAL EXAMPLE 



In this Section, the application of the pair-correlated hyperspherical harmonic (PHH) technique to describe p — d 
scattering for energies above the DBT is discussed. In the PHH approach, the trial w.f. is written as 



LqSq 



(80) 



i=l i=l 

The first term ^Pyi describes the system when the two incident clusters are well separated and it is written in the form 



LS 



where 



(1 - e-T^O^^+i Gl(7?o, god,) ± \Fl{vo, qpd,) 
qodt 



(81) 



(82) 



The factor (1 — e 7di)2L+i j^g^g ^^ggjj^ introduced to regularize the function G at the origin, and 7 is taken as a nonlinear 
variational parameter. Its choice is not a critical problem due to the large flexibility of the trial w.f. in the internal 
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region. The range of values 7 = 0.3 — 0.5 fm"-'^ has been found to be adequate in all the cases studied. For di ^1/7, 
i^Aii) coincides with the first term of the asymptotic form (|73|). 

The second term \l/c of the trial w.f. must describe those configurations of the system where the particles are close 
to each other. For large interparticle separations and energies below the DBT, vE'c goes to zero, whereas for higher 
energies it must reproduce an outgoing three particle state. Each amplitude ipcii) is expanded in terms of the PHH 
basis 

a.T.T K 

where (2)p^>^o jg ^ hyperspherical polynomial, defined for example in Ref. Ua^T.T ,k{p) arc the hyperradial 

functions to be determined by the variational procedure and fa{xi) is a correlation function. In the "uncorrelated" 
hyperspherical harmonic expansion (HH) one has fa{xi) = 1, whereas in the PHH approach it is included in order 
to better take into account those correlations introduced by the repulsion of the potential at short distances. The 
correlation function is calculated by solving a two-body Schroedinger-like equation and it goes to one for large values 



of the interparticle distance Xi. In Eq. (83), the sum over a,T,T is truncated after the inclusion of Nc channels. 
Usually all the channels with £a + < imax, where L,„ax = 6 or 7, are included in the expansion. For each channel, 
the sum over the index K of the hyperspherical polynomial runs from iiT = 0, 1, 2, . . . up to a maximum value Ka.T.T- 
The latter values are chosen in order to achieve the required convergence in the quantities of interest. The total 
number of basis functions is therefore 

M= J2 (Ka.T.T + 1) . (84) 

a,T,T 

To simplify the notation let us label the basis elements with the index fi = [a,T,T, K], and introduce the following 
completely antisymmetric correlated spin-isospin-hyperspherical basis elements 

3 
i=l 

which depends on p through the correlation factor and form a non-orthogonal basis. In terms of 'P^{p, fi) the internal 
part is written as 

The "uncorrelated" basis elements V^{^) are obtained from Eq. ( |85| ) in which all the correlation functions are fa{xi) = 
1. It is important to note that also the elements 'P^{Vl) do not form an orthogonal basis, as has been discussed 
in Ref. [ p6| . In particular, those elements having the same grand-angular quantum number G — ia + La + ZiiT 
but belonging to different channels are not orthogonal among themselves. Moreover, some 'P°(S^) with the same 
GjAcScjT quantum numbers are linearly dependent. In Ref. j3^ these states have been identified and ehminated 
from the expansion. 

In the present case, the basis elements V^{p — > cx),f2) reduce to the uncorrelated ones 'P^{Vl) in the asymptotic 
region Vb(i?). Therefore, it will be useful to combine the correlated basis ( |8^ ) in order to define a new basis with the 
property of being orthonormal for p = oo. This can be readily accomplished by noting that 

"PM'(P) ^ J V^{p, n)^V^,{p, n) ^ Nj^^, + ^+ 0(1/ p') , for p ^ 00 , (87) 
where, in particular, 

j dnvlm^Vl,{n) . (88) 

Let us define a matrix U such that U N^^'^ = V is a diagonal matrix with diagonal elements either 1 or 0. The 
T>f^ — correspond to the states 7^°(ri) linearly dependent on the others. New uncorrelated and correlated bases are 
defined as: 
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M M 

Q° (f]) ^ U^^^'Ki^) ' ^) ^ E ^) . (89) 

^' —1 fi' —1 

The basis functions Qfj,{p, ^) are still not orthogonal for any finite values of p. When p oo, the elements Qfj,{p, ri) — > 
Q°(ri). Due to the fact that some of the uncorrelated elements 7'°(r2) are linearly dependent, some elements Qjl(f^) 
are identically zero. Therefore, some correlated elements have the property: Q^{p, ri) — > as p — *■ oo. 
In terms of the new basis, the internal part '^c is simply 

M 

where = X)^/ (f/"^)^.',x u^>. 

The variation of the functional (|7^) with respect to the hyperradial functions uj^{p), which are the unknown 
quantities entering in the description of the internal part of the w.f. \l/c, leads to the set of inhomogeneous second 
order differential equations given in Ref. [Q, where the method of solution has been discussed. For p ^ oo, neglecting 
terms going to zero faster than the asymptotic expression of such a set of differential equations reduces to the 
form 



where the matrix c is defined as 



+O(J-)^^^,(p) = 0, (91) 



/ dnoPin)^ xQU^) ■ (92) 



In practice, the functions uj^{p) are chosen to be regular at the origin (w^(0) = 0) and, in accordance to Eq. (]9l]), to 
have the following behavior 



M 



M - E (e^'^'"""'),.,' '^^oso,.' e'"^" . (P - oo) , (93) 



where S^^^g^^^, are unknown coefficients. Eq. ( |93| ) reproduces the behavior of the breakup term with 

3 M 



E E ^-.tAP, n,)ya.T.T{^ = E '5loSo,m2m(P, ^) > (94) 
a,T,T 4=1 n—1 



^a,T,r(P:^^0 = E E '^Lo.So.M't^A^V | fc.{x^)^^^Pt''^{e^) , p = [a,T,T,K] , (95) 

K=0 \m' = 1 



and 



M M 

X = EE 2^(^)^ W 2m'(") • (96) 

From the above equation it is clear that an enlargement of the PHH basis by increasing Nc and Ka,T,'T improves 
the trial w.f. both in the internal and in the asymptotic regions. The choice in Eqs. (93) and ( |9^ ) is natural in the 

calculation scheme outlined here and ensures that the matrix element (^'i(,Sol(-ff ^ LoSo) be finite . In fact 

< ^i^k 1^ - ^l*ioS„ r dp I(p) , (97) 



and for p — > oo one has 
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Ap) 



11=1 u' = l ^ ^ ^ 



+c.,(p) c.,,(p) + 1 ^ . (98) 

P P } P 

In conclusion, the procedure appears to be well founded and in the (hypothetical) limit of an infinite expansion it 
leads to the exact results. The breakup amplitudes Aa{p,^i) are obtained in first order approximation in terms of 
the PHH basis elements Qfi{p, i^), via Eq. (p5|). The convergence of A{p, VLi) as a function of Nc and Ka^T,T has been 
found to be very slow when 9i Tr/2. This is related to the difficulties for the expansio n (^5| ) of reproducing the Xi 
dependence of the breakup amplitudes in the region Vi{R) + Vt^i{R) specified in Eqs. (|63|-|65D. However, this problem 
has no consequences on the calculation of the elastic S-matrix elements [S], since the regions [Qi] and [flt.i] do not 
give any contribution to the error term < e\H — E\e >. 

Once proved that the procedure can lead to the correct solution of the problem, there are two main questions to 
be taken care of in practical applications, namely: 

1. The boundary conditions given in Eg. ( ^ ) are valid only at p = oo. However, such a problem can be overcome by 
using the procedure proposed in Ref . Ig^l • First of all, a set of M coupled differential equations are derived from the 
variational condition (5t»jj, [5] = 0. The functions uju,(p) are determined in the region p > po hy expanding in powers 
of 1/p and verifying the boundary conditions of Eq. (p3|). Values of pq w 100 fm have been found to be appropriate. 
Then, the differential equations are solved numerically in the region (0, po) imposing the continuity of the logarithmic 

derivative at p = po with the solutions found for p > po. This procedure fixes the values of the coefhcients Si^^^g^ ^. 
The stability of the solution has been tested by choosing different values of the matching radius in the range 80 — 200 
fm. Finally, the variational condition S-^ [S] = is used to find the elastic coefficients SlqSo ls- 

2. The convergence of the functional [S] with Nc and Ka,T,T- In all the calculations performed pTl, the convergence 
pattern of [S] above the DBT was found to be rather similar to that observed below the threshold To give an 
idea of the convergence, the calculated S-matrix elements for p ~ d scattering at Eiab = {3/2)E — 5 MeV and 10 



MeV are reported in Table VI. The state considered is = 1/2+, which is in general the most structured one since 
it must satisfy the orthogonality condition with respect to the three-nucleon bound state. In this case the S-matrix 
is a 2 X 2 (complex) matrix corresponding to the p — d asymptotic states L — 0, S — 1/2 and L — 2, S — 3/2. The 
NN interaction considered is the AV18 potential ||37|] . 

In Table |Vl| , the real and complex part of the two phase-shifts ^6*1/2, and the mixing parameter ?7i/2+ are 

given for different values of the number of channels Nc = 8, 18, 26. In the first case (Nc = 8) only the 8 channels given 
in Table |l| have been considered since they give the major contribution to the bound state. The cases 18 and 26 
correspond to including all the T = 1/2 channels with ia + La < A and 6, respectively. In all the cases, the number 
of hyperspherical states for each channel has been increased until the convergence was reached. It has been found 
that a number of states similar to the one used for energies below the DBT is sufficient to get converged results, i.e. 
8 basis elements for the first 4 channels, 6 basis elements for channels 5—14 and 4 basis elements for the successive 
channels. 

Also the convergence with Nc is analogous to that obtained for energies below the deuteron breakup, as can be seen 
in Table The first 8 channels give the most important contribution to the S-matrix, with results accurate up to 
three digits. Higher order channels give minor contributions, of the order of a few percent. The contribution of the 
channels with T — 3/2 has been found to be very small and it has been disregarded in the calculation presented here. 



VI. CONCLUSIONS 



The purpose of this paper is to demonstrate the validity of the KVP to describe p — d elastic scattering above 
the DBT. The practical application of the principle has been analyzed, with particular reference to the use of the 
PHH expansion technique. It was shown that the principle remains formally unchanged when the breakup channels 
are open. In this case, the asymptotic behavior of the w.f. has been extensively discussed with special care in the 
treatment of the long range Coulomb potential. Different asymptotic regions has been introduced for different values 
of the hyperspherical variables (p, 9i) and the solution of the Faddeev-Noble equations has been studied in each case. 
It was found that certain asymptotic configurations in the p + p + n system have negligible probability. For example, 
the breakup w.f. vanishes in region Vi when particles j, k are protons or when they are in a relative state with angular 
momentum £a > 0, independently of the fact that they form ap — p or an n — p pair. The transition from the solution 
in region Vb to these vanishing solutions has been studied. 
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Moreover, it was explicitly shown that the use of "approximate" terms in the asymptotic logarithmic distortion 
of the breakup wave function does not influence the formal derivation and the practical application of the principle. 



This is due to the fact that exp(— ixln(2Qp)) x Aa,T,T is in any region a smooth function (as discussed in Sect, pi E| ) . 
Thus it is not very critical to try to reproduce it with some approximate method like an expansion over HH functions. 

The hyperspherical basis has been found to be a rather natural one to describe the asymptotic behavior of the 
wave function describing scattering above the DBT, including also the long-range Coulomb interaction. In fact, the 
equations to be solved have a very simple structure in terms of p in the asymptotic region and, accordingly, analytical 
solutions can be obtained. In practical applications it is not difficult to reduce the problem to the determination of the 
hyperradial functions in the range < p < po with the application of specific boundary conditions at pQ. Appropriate 
values of the matching radius po have been found to be in the range of 80 200 fm. 

The convergence of the quantities of interest is the usual problem which one has to face when applying a variational 
procedure. In two preceding papers pO| , pT[ , it has been shown that the convergence of the elastic observables can 
be achieved with the same precision as below DBT. Realistic NN and 3N interactions can be used to calculate the 
observables in the n ~ d and p — d case, so meaningful comparisons with experimental data can be performed. In the 
present work a detailed discussion of the problem related to possible inaccuracies introduced from an approximate 
description of the outgoing breakup waves was presented, including a numerical example. 

In conclusion, the KVP provide a formalism from which it is possible to obtain reliable and fully-converged results 
for the elastic part of the p — d S-matrix at energies above the DBT even with the inclusion of the long range 
Coulomb interaction. Applications of the method described here were already reported in Rcfs. |^ , ^ where the 
elastic differential cross section and the vector and tensor analyzing powers were studied at Eiab = 5 MeV and 10 
MeV. The further problem related to a precise determination of the breakup amplitudes Aa,T,T is currently under 
progress and it will be discussed elsewhere |Q . 

APPENDIX: 

In this appendix, the asymptotic behavior of the (outgoing) breakup wave of ^(i) in Vi{R) is given. As discussed 



in Sect. IIIB, only the £a — channels of the components '4>j{i) and '(/'*;(*) do not vanish in Vi{R) and we have 



M^) = C{p) i-fl^^^yo^^s^{t) , region V.(i?) , (Al) 



Vfc(*) = C(p)^^^^f^^3^0AASE(*) 



ASS 



where C(p) and ye,L,A,s,^ii) are defined in Eqs. ( p2[ ) and (|T^, respectively. The sum over the channels has been 
restricted to those with £ = as stated previously. The factors {—)^ and l/-\/2 have been introduced for convenience. 
The functions / and g are determined as explained in the following. First of all, (/ase = /ass since ^^ii) must be 
antisymmetric under the exchange j <-> k. Therefore, the amplitude = S.^ Sj + V'fc(*) '^k can be cast 



in the form 



m = C{p)Y.^^^^yAS^i^), (A2) 



A5S 



where 

3^A5s(l) = 5^0AA5s(i)|(t^'t')T,o(t^) 1 +1 ) , T^1~S . (A3) 



The amplitude ip{i) given in Eq. (A2) is easily seen to be antisymmetric with respect to the exchange j <-*■ k. 



The amplitudes and tp{k) in the region Vi{R) can be obtained in the same way as in Sect. [IID, where the 
behavior of in Vj{R) and Vk{R) is discussed. In this region the following relations hold 

Xj ^ -Ky, ^ -iy,; , &j ^ ^ , (A4) 

and 

1 TT 

Xfc ^ KYi Yk -^Yi , - . (A5) 
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Accordingly, ^pj or k{j) and V'j or k{k) satisfy the FNE with y^N _ q ^^^^ coincide with the "free" solution ([24|). On 
the other hand, as discussed in Sect. Ill D , the components V'i(j) and tpi{k) vanish in Vi{R) due to the Coulomb 
repulsion. Therefore, in region Vi{R) 

m = C{p)Y.Al^{^)y^{j) E,+Ai,i^)y^{j) , (A6) 

a 

m = c{p)Y,Al^{^)y^{k) E,+Al,{^)y^{k) , (A7) 

a 

where a = {i, L, A, S, E}. 

Since, under the exchange j <-> fc, it must be verified that t/'O) ~* ~i'{k) and vice versa, the amplitudes 
■^a'j^ m fcC""/^) ^^"^ ^'-'^ independent. Using the previous relations for the Jacobi coordinates, then 

yaU) ^ Y.^~f D^^s'yoAAS'i:{i) , (A8) 

S' 

Mk) ^ ^(-)'+^i?a,s'5^0AAS's(«) , (A9) 

S' 

where 

^-- = H^"''^'^^^'(oo o){M^'}' ^^''^ 



and £ = \/2£ + 1, etc. Note that only states with orbital angular momentum £ — enter the r.h.s of Eqs. (AS) 
and (A9). Finally, we can write 

^(j) + iPik) - C{p) ^AS'S yAS'^i^) , region V,(i?) , (All) 

AS'E 



where 3^as's(«) is given in Eq. (A3) and 

/Cas^s = Efe(?) + ■ (A12) 



ILS 



Inserting the expressions (|Aj) and ( All ) in the FNE and disregarding all terms 0{l/yi) or 0{l/p), the following 
equation is obtained 

IasA^i) + V2S+l{Xi)fASY,iXi) = -X,V2S+liXi) /Cass , (A13) 

m 

where 

dx, dy, (3^A'S'S'(«)) VnN(x,)3;ass(«) = V2s+i(a;0<5ss"^AA'feE' ■ (A14) 

V2s+iixi) is the projection of the NN potential acting in a state with orbital angular momentum zero, spin S and 
isospin T — 1 — S. 

The general solution of Eq. ( A13) can be written as the sum of the solution of the homogeneous part and a particular 
solution of the complete equation. The solution of the homogeneous part will be denoted as jAST,f2S+ii^i)' where 
7ASi; is an arbitrary coefficient; f2s+i coincides with the (regular) solution of the zero-energy Schroedinger equation 
of a two-nucleon system in the 25+1 s-wave state. For large Xi values it goes like Xi — 0,23+1, where ai (03) is the 
singlet (triplet) scattering length. 

A particular solution is simply 

/ass = ^^SY, ■ (A15) 
Note that the particular solution exactly cancels the amplitudes + tp{k) in Vi{H). In fact, 

* = tp (i) + tPi j) + ij{k) , 

« C{p) E 7A5s^%^3^A5e(*) . (A16) 

ASS ' 
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Namely, in Vi{R) the pair jk is in a relative s-wave state with vanishing relative velocity, whereas the motion of the 

third particle i is governed by the distorted wave C{p) « C{yi). 

For large values of Xi, fhST,{xi) / Xi — > const. In the region Vt.i{R), the amplitude given in Eq. ( [Al[ ) should smoothly 
match the behavior of ipk{i)- Therefore, the constants 7ass are related to the breakup amplitudes A^^ j via 

the relations 

^(7ASE - /Cass) = IU5E,,('r/2) . (A17) 
Finally, the functions (j^ai^i) introduced in Eq. (E^) have been defined to be 



^ ^?)AASE,,(^/2) 



''OAASS 
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TP, ^ — ^ A/TpV 




N, = 8 


M — 1 8 


/V — 9fi 


"Dm 




(-5.45,0.004) 


(-5.43,0.004) 


(-5.43,0.004) 


Sl/2 




(-42.0,1.74) 


(-41.8,1.74) 


(-41.8,1.74) 


V1/2+ 




(1.06,-0.03) 


(1.06,-0.03) 


(1.05,-0.03) 


Eiab = 10 MeV 




iVc = 8 


= 18 


A^; = 26 






(-7.41,0.23) 


(-7.31,0.24) 


(-7.30,0.24) 


Sl/2 




(-60.8,11.7) 


(-60.6,11.7) 


(-60.6,11.7) 


V1/2+ 




(1.04,0.06) 


(1.02, 0.06) 


(1.01,0.06) 


TABLE I. Convergence 


of the oigcnphaso shift and 


mixing angle parameters of elastic p — 


d scattering as a function of the 


number of channels included in the PHH expansion of the internal part of the w.f. The state considered has quantum number 


— 1/2^. The nuclear potential model considered 


is the AVIS interaction. Eiab is the laboratory energy of the incident 


proton (the center of mass 


energy is E 


= 2/3Elab)- 






a £ct 


Let 


Ac 


Sc T 


E„ T 


1 








1 


1/ 2 1/2 


2 








1 


1/ 2 1/2 


3 


2 


2 


1 


3/ 2 1/2 


4 2 





2 


1 


3/ 2 1/2 


5 2 


2 





1 


1/ 2 1/2 


6 2 


2 


2 


1 


3/ 2 1/2 


7 2 


2 


1 


1 


1/ 2 1/2 


8 2 





1 


1 


.'!/ 2 1/2 



TABLE II. Quantum numbers for the first 8 channels considered in the expansion of the internal part of the three-nucleon 
scattering w.f. 
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